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Spin critical opalescence in zero temperature Bose-Einstein Con- 
densates 
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Abstract -Cold atom developments suggest the prospect of measuring scaling properties and long- 
range fluctuations of continuous phase transitions at zero-temperature. We discuss the conditions 
for characterizing the phase separation of Bose-Einstein condensates of boson atoms in two distinct 
hyperfine spin states. The mean-field description breaks down as the system approaches the 
transition from the miscible side. An effective spin description clarifies the ferromagnetic nature 
of the transition. We show that a difference in the scattering lengths for the bosons in the same 
spin state leads to an effective internal magnetic field. The conditions at which the internal 
magnetic field vanishes (i.e., equal values of the like-boson scattering lengths) is a special point. 
We show that the long range density fluctuations are suppressed near that point while the effective 
spin exhibits the long-range fluctuations that characterize critical points. The zero-temperature 
system exhibits critical opalescence with respect to long wavelength waves of impurity atoms that 
interact with the bosons in a spin-dependent manner. 
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Introduction. — Now that cold atom technology has 
realized uniform trapping potentials bounded by sharp 
edges [IJ H], the scaling of near zero-temperature phase 
transitions can be explored in the laboratory. We con- 
sider spin domain formation in a Bose-Einstein condensate 
(BEC) of atoms in two distinct hyperfine states, which 
we call the "spin up" (|f)) and "spin down" (\\)) states. 
The Feshbach tuning of one of the scattering lengths can 
trigger this transition in the quantum (zero-temperature) 
regime. In-situ images of the atoms in one of the spin 
states can reveal spin density fluctuations at a character- 
istic length scale (the correlation length) that diverges as 
the scattering length is tuned near a critical value. We 
show that the mean-field description breaks down near 
the transition so that the critical exponents may differ 
from their mean-field values. To provide a reference to 
future experiments and to reveal trends with respect to 
polarization, scattering lengths, and density, we calculate 
fluctuation properties in mean-field. Using a spin descrip- 
tion, we find that, for Ising spin-spin interactions with 
equal j" — j" and I — l scattering lengths, the spins ex- 
hibit long-range fluctuations whereas the long-range den- 
sity fluctuations are suppressed near the transition. The 



system then remains transparent to distinguishable low 
energy atoms that interact with the bosons in a spin in- 
dependent manner whereas it turns opaque to atoms that 
experience a spin-dependent interaction. This system can 
be realized with 87 Rb atoms that support a resonance in 
the interaction of atoms in different spin states [3J|3]. 

Switching ground states. — The phase separation 
of boson supcrfluids was predicted in [5], its cold atom 
realization, dynamics and surface tension were predicted 
and described in [5] and the transition was demonstrated 
in [7J. Here, we show that the ground state configuration 
of a dilute homogeneous BEC of N± spin-down bosons 
and Nf spin-up bosons of mass m confined to a macro- 
scopic volume f2 alters from a homogeneous 'miscible' mix- 
ture to an 'immiscible' separated state. Bosons at posi- 
tions x and x' interact via effective short-range potentials 
Aj.(-|-)5 (x — x') if they occupy the down (up) spin state 
and via XjjS (x — x') if their spins differ. The interaction 
strengths A arc proportional to the respective scattering 
lengths, \+ iu — {^irh 2 /m)a^^u , one of which can be 
Feshbach-tuned across the phase boundary. 

In the immiscible ground state configuration, 7Vj_ spin- 
down bosons reside in a volume f2j, and Nf spin-up bosons 
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occupy fif = Q — Assuming the volumes are suf- 
ficiently large to neglect surface effects, the immiscible 
ground state energy is 



E 



sep 

We determine Oj, 



2n 



A 



2\n 



(i) 



by minimizing E sep , which gives 
XlNf/ 2Vt\ = A t iV|/ ^(O-Qj.) 2 ], implying equal 
pressures in £1 i and O-j-. Inserting the corresponding vol- 
ume fraction, fl/fli = I + (iV-f /-/V4.) a/ / Aj, into E'sep, 



ft 



20j 



•(0-%) 



A t iV T 2 



A, A^p 
= 04—1, ( 2 ) 



we find that the separated ground state energy 
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drops below the mean-field energy E m i X of the homoge- 
neous mixture, 



E,, 



^ 2 + ^ + ^», 



(3) 



if yAfA^ < Ay. Characterizing the competition between 
like and unlike interactions by 



AfAj, 

the BEC ground state switches at <? = 1. 



(4) 



Divergence of single component compressibility. 

— We determine the spin-down density response Spi to 
a potential perturbation SV^. The low frequency, long 
wavelength response follows from the Thomas-Fermi de- 
scription that minimizes the homogeneous free energy 
Fq = E mix — fJ-iNi — fi-^Nf at fixed chemical potentials 
Hi and and replaces pi ~ ► P-i ~ <5^V We solve the 
resulting equations 



hPi + x uPt = AH - SVi 



(5) 



by linearizing Puj-) — Pitf),o + $Put) around the homo- 
geneous (SVi = 0) equilibrium densities Pim,o- The 5pi 
result gives a response function (proportional to the com- 
pressibility) 



SVi 



h [1 - 9] 



(6) 



that diverges when g — > 1. The divergence of the com- 
pressibility implies large scale fluctuations that cause crit- 
ical opalescence [5] near ordinary critical points. The con- 
tribution of [— (Xig)} to the x|-denominator of eq. (|6|) de- 
scribes the long wavelength 4.-4. attraction mediated by the 



t" BEC that competes with the short-range 4--4- repulsion 
described by Aj.. The divergence implies that the small- 
est of Vj.-potential variations induces a large pi density 
response. Equating the flat potential requirement for sim- 
ulating infinite systems to Spi/pi^ <C 1 suggests that Vj, 
control should ensure that 

AV^«^[l- 5 ], (7) 

where AV4. represents the 4- potential variation over fl. 

Mean-field breakdown. — The mean-field descrip- 
tion predicts the boundary of its validity regime [9] . The 
misciblc (g < 1) mean-field densities are homogeneous, 
except for small amplitude fluctuations. This assump- 
tion breaks down when (^-quantum fluctuations induce 
an effective SVi — At/Jp-^-variation sufficiently large to 
give \5pi\ ~ Pi.o- Estimating the quantum fluctuation 
as \Spi\ ~ ^/lim r 

1/4 



+o(6pf (r) Sp t (0)) so that \5p t \/p ^ 
we expect large Spi fluctuations when 



I 3\V4 



(8) 



where a-j- (aj,) is the scattering length of spin- up (down) 
state. A similar relationship follows for Sp^. Assuming 
similar densities and scattering lengths, p^ ~ pi ~ p/2, 
where p is the total density, p = p^ + pi, a-j- ~ aj. ~ ajj 



1 /3 

and introducing a = (ofa^aj/) we suggest that 



[1 - 5] > (pa 3 



,V4 



(9) 



is a necessary condition for mean-field validity. 



Feshbach steering of the system and required 
magnetic-field control. — In a magnetically controlled 
Feshbach resonance it is the strength of a homogeneous 
magnetic field B near the resonant value Bo that varies 
the scattering length. How accurate does the magnetic 
field have to be controlled to avoid the effects of a fluctu- 
ating interaction? The magnetic field dependence of the 
scattering length 



a = ciq 



1 



B-B 



(10) 



where ao denotes the background scattering length and A 
represents the width, implies a scattering length variation, 
Sa, induced by a magnetic field variation SB equal to 



5a = ao 



A 



B-Bn 



SB 
A"' 



(11) 



For the S7 Rb case, the near-equality of the triplet and 
singlet scattering lengths leads to a-j- ~ aj, ~ ai ~ 5nm 
within a few percent. Hence B should be tuned far from 
resonance, \B — Bq\ ~ A/S to encounter the phase bound- 
ary. This is important as the resonance of refs. [3], [4] has 
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been found to quite 'lossy'. At \B - B \ ~ 50A - 100A, 
however, particle loss should not play a role. This im- 
plies that the transition condition is achieved far from the 
resonance with B A B() — S a few percent. The relative 
scattering length variation is, then, of order <5 2 even if 

|<ys|~A, 



Sa 


= S 2 


SB 


a 




A 



(12) 



Hence, it should be feasible to ensure that \Sg/g\ = 
2\Sa u /a u \ « [1-g]. 

Spin analogy. — Scaling exponents have been de- 
termined for finite temperature transitions of spin lat- 
tices. We introduce the effective spin operator a so that 
o z ||) = + ||), a z ||) = — ||), and the a x and <r y -operators 
are represented by the Pauli-matriccs in the ||), ||) basis. 
The interaction of bosons with coordinates r and r', where 
r and r' represent both location x and x' and spin a and 
is [U3] 

V eff (r, r') = X ° [(A t + A, + 2X V ) 

+ (A t - h) + K) 

+ (A t + X l - 2A„) a z a' z \ . (13) 

The term linear in a z describes an effective short-range 
magnetic field carried by the particles. This effective field 
interacts with the other spins. The corresponding inter- 
action term contributes a mean-field energy that is indis- 
tinguishable from that of an effective magnetic field h e ff, 
h-eff ■ o") where h e ff = [(At — Aj,) /4]pz. Characterizing 
the interactions by 

Af + Aj. , Af - Xi Xu , 1<n 

Al = ^^' d= MTTi ' and r= A7' (14) 

the spin-spin interaction potential reads 

K//(r,r') = <5(x-x')^[(l + r) 

+ d(a z + a z ) + (l-r)(a z cr z )}, (15) 

g takes the form g = r 2 /[l — d 2 ] and h. e ff = {d/2)Xi,pz. 
The zero internal magnetic field condition, d = 0, provides 
a special point: only for h. e ff = do the ground state 
spins locally align when the spin-spin coupling term turns 
ferromagnetic as illustrated in fig. [TJ The spin analogy 
suggests that (cr z ) plays the role of the order parametciQ 

1 Measurements of (& z ) in the f2-j- and Q^-volumes would record 
the different branches of the magnetization curve. At finite tem- 
perature, {<j z ) could vary discontinuously across the transition as a 
function of T, except for P = 0atd = 0- the critical point. In anal- 
ogy with finite temperature phase separation transitions, we suggest 
that there may be a line (plotted as a function of P) of spinodal de- 
composition and a coexistence line that touch at the critical point. 
The second-order nature of the zero-temperature transition (inde- 
pendent of P) may be a consequence of these lines approaching each 
other as T — > 0. At P = 0, d = 0, and fixed temperature value, the 
(& z ) -variation as a function of r tends to its asymptotic value ±1 
over an r-interval of magnitude kgT/ [X^p], so that for T = 0, (a z ) 
'jumps' from to ±1, but this quantity still provides a legitimate 
order parameter. 




Fig. 1: In spin-language, the difference of || and || scatter- 
ing lengths translates into a short-range internal magnetic field 
carried by the particles interacting with the other spins. This 
internal magnetic field can locally align the effective ground 
state spins, even if the spin-spin coupling is antiferromagnetic 
(Xl (1 — r) < 0). In this figure, the at, aj. phase diagram for 
fixed a u value is shown. We illustrate the — effect on 
the zero-temperature transition by showing that the immisci- 
ble regime (shaded by horizontal lines) and the ferromagnetic 
coupling regime (shaded by diagonal lines) only intersect along 
the at = aj.-line. Only when au is varied while a-f = aj., does 
phase separation at zero temperature take place when the spin 
coupling switches from ferro- to antiferromagnetic. 

possibly offset by its value in the homogeneous mixture. 

Zero temperature itinerant ferromagnet-like 
transition while avoiding "lower branch" physics. 

— The effective spin description reveals that the tran- 
sition combines ingredients of phase separation and fer- 
romagnetic transitions |11) . The fermion analogue, dis- 
cussed in the pioneering work on quantum phase transi- 
tions [12], was reported in a cold atom trap [T3]. That 
transition, however, requires strong repulsion as the com- 
ponents separate when the inter-particle interaction en- 
ergy outweighs the kinetic energy. In this strong coupling 
regime, two-particle bound states (dimers) form involv- 
ing 'lower branch' physics P3HH]. In contrast, BEC- 
phasc separation is triggered by the competition of dif- 
ferent short-range interactions, all of which can be weak, 
allowing the system to remain in its metastable BEC state. 

Mean-field fluctuations. — We describe the correla- 
tions induced by the quantum fluctuations of the normal 
modes that diagonalize the free energy operator 

F = H — fH N t - mN h (16) 

where H denotes the Hamiltonian. The zero-momentum 
replacement of the creation and annihilation operators 
Cj,k=u ; Cj k=0 — > y/Nj where j =|,| leads to the above 
mean-field free energy 

F Q = E mix - fi t N t - mNi. (17) 
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The next order contribution in the Bogoliubov expanded 
free energy F = SF + Fq keeps the quadratic terms 



SF= J2 ( S J,k 6 i>k 




where = h 2 k 2 /2m, and the chemical potentials, fj^-, 
Hi, were replaced by eq. ([5} with SVi = 0. Two- 
component BEC Bogoliubov transformations were con- 
structed long ago, but these treatments obscure the un- 
derlying oscillator-structure [20], [21]. As in Ref. [10] we 
introduce the phase space position/momentum-like oper- 
ators 



Vj,k 



b ],k + h-* 
V2 



t\ j^t A (19) 



that map SF onto a sum of oscillator Hamiltonians. Be- 
low, the fc— subscript and the sum over repeated indices 
will be tacitly understood. Unlike Ref. [10], we make use 
of the simplectic approach to phase space transformations 



[23] . We denote the phase space vector ( 
and write SF = (l/2)QKijQ where 



K 



K, 









i^,1Tl) 

(20) 
(21) 

The Bogoliubov transformation is a linear, canonical point 
transformation £ — > rj, rj = (</>+, <fi—, 7r+, 7r_) and Q = 
MijTjj where 



K _ ( e + A t n t K^/nyni 
Xus/n^nl e + A t n t 



,K n = 



M 



M 



M, 

resulting in a transformed K-matrix of the form 



K' = MKM 



E 
E 



(22) 



(23) 



where M represents the transpose of M and E is the di- 
agonal matrix of collective mode eigenvalues 



E 



E+ 



E_ 

The eigenvalue equations (f22]) and (f23j) take the form 
E =M,KA, 



(24) 



E 



(25) 
(26) 



The K' is then equivalent (up to a constant) to the di- 
agonalizcd Bogoliubov Hamiltonian SF' = J2a=± ^^a^a 
where b\ = + iir a ) /y/2, b a = {<t> a - in„)/\/2. To en- 
sure that b a ,b^ a satisfy the boson commutator relations, 
the M transformation has to be canonical [23| implying 



(27) 



Writing as the product of a rotation matrix R (R = 
R _1 ) and a diagonal scaling matrix S, we have = RS 
with 



i cos 9 — sin 9 , 
sin 9 cos 9 





r_ 



(28) 



As 



We satisfy Eq. (J27J) by choosing = RS 

Kjr = el, Eq. ([MJ leads to S~ 2 = E/e or T± = y/e/E±. 
Inserting the corresponding S into Eq. (|25p results in 



RK„iR = S : ES 1 = 



E 2 + /e 
E^/e 



(29) 



so that E follows from the diagonalization of and E+ / e 
and E^/e are the eigenvalues of K^. For a mixture of 
overall polarization, P = (JVf — iVj,), we cast the result- 
ing transformation in terms of the average sound velocity 
c = y Xlp/itl, the length £ = K/ (mc) and a transition 
parameter 



t = (1 - .9) 



(1 - d 2 ) (1 - P 2 



(1 + dP) 2 
Diagonalizing the K^, we find 

El =h 2 k 2 c 2 ± (l + e±k 2 ) 



(30) 



where 



_ « 



dp 



(31) 

(32) 
(33) 

with S (x) = and the rotation angle 9 in R is 

determined by 



vT+dP V l±v / T :r t' 
The S matrix elements, r±, take the form F± 



cos (9) 



d + P 



1 



1 + dP J y/T=t 



(34) 



The corresponding transformation fifd) ^+(-)' 

7r ta) n +( _ } with (nt fciv,*,) = (s^iv,^) = 

(l/2)S cr ^ (T 'Sk,k' determines all correlation and response 
functions. For instance, the mean-field tt"density correla- 
tion function is 



(fa ( x ) fa (0)) = Po,t + #>,t 



(2^ 



■ik-x 



(35) 
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where 



2(0 k)t 0_ k , t ) = s (S+k) cos 2 (9) + s(f_fc) sin 2 (0) . (36) 



Note that the correlation functions harbor two length 
scales: £ + and As diverges, £_ sa 



i~ 1/,2 £y 2/ (1 + dP) 2 , this length scale, the correlation 
length, exceeds the imaging resolution before the tran- 
sition is reached. 

For the special case of Ising spin-spin interactions, d = 
0, we determine the long-range part of the density-density 
and spin-spin correlation functions, 



(P(x)p(0)) ~ P l + p a 



(a z (x) a z (0)) 



d 3 k 
(2^7 



e^'sit-QA- (37) 



Po 



d 3 k 



(38) 

where po = Po.t + Po,i an d where the ^4 -amplitudes de- 
pend on the Bogoliubov angle, giving 



p 



s 2 



VT+PWI- 



p 



(39) 



.4 



p 



v 2 



vi + pwi- 



p 



(40) 



Hence, the near-transition spin-spin correlation function 
exhibits the long-range order parameter correlations typ- 
ical of a second-order phase transition. In contrast, the 
long-range part of the density-density correlations are sup- 
pressed, A 



/>/> 



t if P =/= near the transition. Even 
though the mean-field approximation breaks down, we 
suggest that the near-transition suppression of long-range 
density fluctuations is an actual feature of d = fluctua- 
tions. 

The correlation functions can be extracted from the 
pixel counts of a single in-situ f-density image, averag- 
ing pf (R + x/2) pf (R — x/2) over R (making the self- 
averaging assumption that the R-average is equivalent to 
averaging over many samples). In accordance with the 
suppression of long-range density fluctuations, we suggest 
that for a single image can be extracted from the 'f- 
density. A p^ image that resolves features on the £_ length 
scale though not on the scale can be converted into a 
Pl image by pi ps po — pf. In the case of 87 Kb, d nearly 
vanishes but is not exactly zero (d ~ 0.01). Therefore, 




-0.03 -0.02 



-0.01 0.00 0.01 



0.02 0.03 



Fig. 2: Plots of the ratio of the long-range amplitudes of 
the zero-temperature density-density and spin-spin correlation 
functions as a function of the relative difference d of the ft an d 
scattering lengths, d = [a-j- — aj / [a-j- + Oj.], for fixed coher- 
ence length £_ = 10£ and different values of the polarization, 
P = [Nf- Ni] /N, with P = (full line), P = 0.2 (dashed 
line), P = 0.5 (dotted line) and P = 0.7 (dash-dotted line). 



we investigate the effect of a small but finite d-value on 
the ratio of the long-range density-density and spin-spin 
correlation amplitudes in fig. [2j In this figure, we plot 
the A~ p /A~^ az ratio as a function of d for different P- 
polarizations when £_ = 10£. Note that for d ranging 
from —0.03 to 0.03, the A~ p /A~ r(J __ -ratio remains smaller 
than one part in one thousand with P ranging up to 0.7, 
suggesting that the near transition suppression of the long- 
range density fluctuations remains valid even for small but 
finite d- values. 

Spin opalescence. — In a classical system in equilib- 
rium at temperature T, the long wavelength structure fac- 
tor is proportional to koTn where k denotes the isother- 
mal compressibility. Near the critical point of a gas-liquid 
transition, the paradigm of critical opalescence, k diverges. 
The compressibility is also the long wavelength limit of 
the static density response function. In the ground state 
of the spinl/2 BEC system with Ising spin interactions 
(d=0), it is the spin response that diverges, not the den- 
sity response. To show that, we consider a weak mag- 
netic field perturbation of good momentum k described 
by a contribution (5-ffe lk ' x £ z (x)/f2 to the energy density 

where £ z (x) = ^(x^^x) — ^(x^^x) /Q to the en- 
ergy density. The long-range part of the mean-field static 



response <5£ z (x) 
where 



SE z e lk - x /fl is given by SE Z 



\ r. 



Xa-a, = A o 



1 + k 2 C 



SH 



(41) 



which diverges cx t^ 1 at the transition, whereas the 
analoguous long-range density response function Xp P 
A~ p /c?_ cx t vanishes near the transition. As a conse- 
quence, the spin fluctuations can mediate interactions in 
a very pronounced manner and the Hamiltonian terms ne- 
glected in the Bogoliubov approximation become impor- 
tant. 
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Feasibility and summary. — Critical slowing pre- 
vents an actual zero-temperature crossing of a second- 
order phase transition: the temperature T should be lower 
than an energy scale, mc 2 _, that vanishes at the transi- 
tion (mc 2 _ ~ t in mean- field). Experimentally, it would be 
very interesting to reach the mean-field breakdown regime 
where the fluctuations acquire large amplitudes, under the 
conditions of experimental access: ksT < mc 2 _, the mea- 
suring time is tm S> h/[mc_], the potential variation AV 
remains sufficiently small, and the system's linear size L 
significantly exceeds the coherence length. For a 87 Kb two- 
component BEC of density p ~ 5 x 10 13 cm~ 3 , P ~ 0.1, 
a ~ 5nm, we estimate that mean-field breaks down at 
g ~ 0.95, where £ ~ 5£. The above values give L ^> 5£, 
k B T < mc 2 _/20, t m > 20ft/[mc 2 ] and AV < mc 2 /20. 
These conditions, while challenging, can be met in exper- 
iments. 

We have discussed the prospect of characterizing the 
dilute gas BEC phase separation of bosons in distinct hy- 
pcrfine states as a zero-temperature second-order phase 
transition. For equal like-boson scattering lengths, we ex- 
pect the system to exhibit long-range spin-spin correla- 
tions whereas the long-range density fluctuations are sup- 
pressed. 
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